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ABSTRACT 

We spectroscopically study the variability of Ha emission in mid- to late-M 
dwarfs on timescales of ~ 0.1 — 1 hr as a proxy for magnetic variability. About 
80% of our sample exhibits statistically significant variability on the full range 
of timescales probed by the observations, and with amplitude ratios in the range 
of ~ 1.2 — 4. No events with an order of magnitude increase in Ha luminosity 
were detected, indicating that their rate is < 0.05 hr^^ (95% confidence level). 
We find a clear increase in variability with later spectral type, despite an overall 
decrease in Ha "activity" (i.e., Lna/i^boi)- For the ensemble of Ha variability 
events, we find a nearly order of magnitude increase in the number of events 
from timescales of about 10 to 30 min, followed by a roughly uniform distribution 
at longer durations. The event amplitudes follow an exponential distribution 
with a characteristic scale of Max(EW)/Min(EW) — 1 0.7. This distribution 
predicts a low rate of ~ 10~^ hr~^ for events with Max(EW)/Min(EW) > 10, but 
serendipitous detections of such events in the past suggests that they represent 
a different distribution. Finally, we find a possible decline in the amplitude of 
events with durations of > 0.5 hr, which may point to a typical energy release in 
Ha events for each spectral type (-Ehq ~ x t ~ const). Longer observations 
of individual active objects are required to further investigate this possibility. 
Similarly, a larger sample may shed light on whether Ha variability correlates 
with properties such as age or rotation velocity. 
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Introduction 



One of the primary indicators of magnetic heating and activity in low mass stars is Ha 
chromospheric emission, which traces the presence of gas at temperatures of ~ 5000 — 10, 000 
K. In M dwarfs and later spectral types Ha is a particularly prominent indicator since it is 
more easily accessible than other chromospheric li nes such as Ca II a nd Mg II, which are 
located in the faint blue part of the spectrum (e.g.. iHawley et al.lll996l ). Moreover, informa- 
tion on the Ha line is readily available as a by-product of any spectroscopic observations that 
are used to classify M dwarf properties, such as the TiO band-head at 7050 A . Thus, large 
samples of M dwarfs now exist with measurements of Ha luminosity and its ratio relative to 
the bolometric luminosity, LnaZ-^^boi (commonly referred to as the Ha "activity"). 

These samples have led to several important results concerning chromospheric activity 
in low mass stars. First, the fraction of objects that exhibit Ha in emission increases rapidly 
from ~ 5% in the K5-M3 dwarfs to a peak of ~ 70% around spectral type M7 , followed 



by a subsequent decline in the L dwarfs (iGizis et al.ll2000l : IWest et al.l 12004 |2008| ). Second, 
while the level of activity increases with both rotation and youth in F-K stars, it reaches 
a saturated value of Lna/i^boi ~ 10 ~^'^ in M0-M6 dwarfs, followed by a rapid decline t o 
Lna/Lboi ~ 10"^ by spectral type LO JHawlev et al.lll996l : [Giz"is et allboool : IWest et al.ll2004[ ). 
Third, a rotation-activity relation is observed in spectral types earlier than ~ M7, such that 
essentially all objects with f sin(i) > 5 km s^^ exhibit saturated Ha activity. However, 
late- M and L dwarfs exhibit reduced activity even at high rotation rates, f sin(z) > 10 km 



s"^ flMohantv fc Basril l2003l : iReiners fc Basril l2008f ). Finally, a small fraction (< 5%) of 



late-M and L dwarfs have been serendipitously observed to exhibit Ha flare s that reach the 
saturated emission levels found in the earlier-M dwarfs (iLiebert et al.ll2003l ). 



In this paper we focus on the last point (Ha variability) in the spectral type range 
M3.5-M8.5, which encompasses the peak of the Ha active fraction, the regime of saturated 
emission, and the breakdown of the rotation-activity relation. In this spectr al type range 
stars are also fully convective, indicating that the solar-type aQ dynamo (e.g.. |Parkerlll955l ) 
no longer operates. Thus, studies of Ha temporal variability provide additional constraints 
on the magnetic dynamo mechanism. While a large number of objects in the mid- and 
late-M spectral t ype range have b een included in the studies outlined above (for example, 
~ 10^ objects in lWest et al.ll2008l ). Ha variability was not a key aspect of the observations 
and thus most of the flare and variability detections have been serendipitous in nature. As a 
result, no systematic results on Ha variability timescales and amplitudes are available in the 
literature from controlled and uniform cadence observations of a large sample of mid- and 
late-M dwarfs. Indeed, the few existing st udies have only targete d Ha variability in small 
samples of early M dwarfs. For example, iBopp fc Schmitzl (119781 ) studied 15 objects with 
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spectral types earlier than M4.5 and with a cadence of > 1 day, whil e Pettersen et al.l (119841 ) 
studied Ha variability in only three M3-M3.5 dwarfs. iGizis et al.l (120021 ) studied a larger 
sample of mid- and late-M dwarfs, but did not use a uniform cadence (or discuss what their 
cadence was). Finally, observations of a small number of late-M and early-L dwarfs with 
durations of ~ 8 — 10 hr have revealed periodic Ha emission in two objects, tied to their 
rotation period (IBerger et al.ll2008l . l2009l ). 

Here we present spectroscopic observations of over 40 M3.5-M8.5 dwarfs designed to 
probe chromospheric variability on timescales of about 5 min to 1 hr. Focusing on the 
objects that exhibit Ha emission, we find that nearly 80% are variable over the full range 
of timescales probed by our observations. The outline of the paper is as follows. In ^ we 
describe the observations and measurements of the Ha equivalent widths and fluxes. We 
study the Ha emission and its variability in §31 and finally discuss the observed trends and 
their implications in §4l 



Observations and Data Reduction 



We targe ted 43 well-studied M dwarfs in th e spectral range M3 . 5 to M8.5, selected from 
the sa mples of belfosse et al.Nl998h.lGizislJ2002h.lMohanty fc BasrilJ2003h.|Phan-Bao fc Bessell 



kood l. M^antv fc Basrrj2003h . ICruz et al.l J2003h . and lCrifo et al.l J2005h . Given the ob- 
servational setup (see below), we selected targets at a distance of < 25 pc and with an 
observed magnitude of < 20 mag. The properties of our selected targets, including lumi- 
nosities and rotation velocities when available, are summarized in Table [H Robust ages are 
not available for our objects, but the small nearby volume of the sample indicates that few 
objects are expected to be very young. The majority of the objects in our sample have been 
previously shown to have Ha emission, although only VB 10 was known to exhibit flaring in 
the chromospheric emission lines. 

The observations were carried out using the Boiler & Chivens Spectrograph mounted on 
the du Pont 2.5-m telescope at Las Campanas Observatory, Chile, on two separate occasions: 
2007 March 14-17 and 2007 September 12-17. In all observations we used a slit width of 
1.5", matched to the average seeing conditions, and a 600 lines mm~^ grating blazed at 5000 
A . The spectral coverage extended from about 3680 to 6850 A , designed to cover a wide 
range of the hydrogen Balmer lines, the Ca II H&K doublet, and the He I lines. The spectral 
resolution was about 5 A . 

The March 2007 observations were carried out in good weather conditions, with typical 
seeing of about 1.2", while the conditions during the September 2007 observations were 
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poorer, with strong winds and a typical seeing of about 1.5". We observed a total of 20 and 
23 objects in the two runs, respectively, with individual exposures of 300 or 600 s (depending 
on the brightness of the object) and a total observing time of about 1 hr per source. Six of 
our objects were observed more than once. In total, we obtained over 600 individual spectra 
spanning ~ 3000 min of total exposure time. A log of the observations is provided in Table [1] 
The data were reduced using standard routines in IRAF, and the wavelength calibration was 
performed using He-Ar arc lamp exposures. 

To measure the Ha equivalent width (EW) in a uniform manner we fit a second-order 
polynomial to the pseudo-continuum from 6540 to 6620 A (excluding ±10 A around the Ha 
emission line). The equivalent widths were then determined by summing the area under 
the Ha line. The error on each equivalent width measurement includes the noise in the 
spectrum over the same spectral range, as well as the uncertainty in the continuum fit (using 
a statistic). 

3. Ha Equivalent Width Light Curves and Overall Variability 

The Ha equivalent width light curves for all individual observations are shown in Fig- 
ures [1] and O We highlight several individual sources in Figure [3l In Table [2] we provide for 
each source a summary of the median, minimum, maximum, and root-mean-square (RMS) 
scatter of the Ha equivalent width (EW). 

We determine whether a source is variable using a test. Namely, we fit a straight line 
at constant EW through the light curve and calculate the value for the best fit. Nine of 
the objects have light curves that are consistent with non- varying Ha emission on timescales 
of ~ 5 — 60 min at a confidence level greater than 95% (one of them, 2M2226— 7503, has 
two separate observations consistent with non- varying emission). These objects are denoted 
with a '(C)' in the column of RMS values in Table [21 Figure [3^ shows the light curve of one 
such object. We thus conclude that only ~ 20% of mid- to late-M dwarfs with Ha emission 
are non-variable on a timescale of < 1 hr. 

The rest of the sources (~ 80%) exhibit a wide range of variability timescales and am- 
plitudes, from rapid variations at the minimum time resolution of our observations (e.g., the 
M5.5 object 2M0253— 7959, Figure (Hb) to slower variations that span the entire observation 
(e.g., the M7 object 2M1309— 2330, Figure [SJ:;). Similarly, the variability amplitudes range 
from about 1 A to over 20 A . Naturally, the sensitivity to small amplitude variations is a 
function of the signal-to-noise ratio, which in turn depends on the brightness and hence spec- 
tral type. The typical variability amplitudes corresponds to fluctuations of about a factor of 
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two in the Ha luminosity. 



We use several indicators to quantify the variability strength. The simplest quantity 
is A(EW) = Max(EW) — Min(EW), the difference between the maximum and minimum 
EW values observed for each object during the course of our observations. However, since 
the conversion between equivalent width and luminosity depends on the pseudo-continuum 
brightness (and hence spectral type), this quantity cannot be easily compared across spectral 
types. To account for the variation in continuum luminosity between spectral types, we use 
the ratio of the maximum and minimum EW values, i?(EW) = Max(EW)/Min(EW), and 
the RMS normalized by the median equivalent width, RMS(EW)/(EW). These quantities 
are plotted as a function of spectral type in Figure |H In all three cases we find a clear 
rising trend in variability as a function of spectral type, with an apparent flattening beyond 
spectral type M7. 

The non-varying objects are included in these and subsequent plots for completeness. 
However, we note that while they may appear to exhibit significant variability as measured 
by these metrics, their error bars are correspondingly larger. In addition, for several objects 
we have more than one observation (Tabled]). We treat the data for these objects as separate 
observations in order to maintain a uniform cadence across our sample, although in Table [2] 
we use the combined data. 

The distribution of RMS(EW)/(EW) as a function of (EW) is shown in Figure O We 
find no clear correlation between the variability and mean equivalent width. The typical value 
of RMS(EW)/(EW) is ^ 0.25, and only about 10% of the objects exceed RMS(EW)/(EW) ^ 
0.5. 

Finally, we investigate the variability in terms of Ha luminosity. The conversion be- 
tween EW and log(LHQ:/i^boi) is a function of the spectral type since the continuum lumi- 
nosity declines with later spectral ty p es. H ere we adopt the conversion values (so-called x 
factor values) from IWalkowicz et al.l (120041 ). The resulting mean and maximum values of 
log(LHa/^boi) are listed in Table [2j In Figure [6] we plot the range of maximum and min- 
imum \og{LYia/ L\,o\) values for each object as a function of spectral type. We recover the 
same overall declini ng trend in the mean Ha activity as a function of spectral type demon- 



strat ed previously (IHawlev et all I2OOOI : ICruz k ReidI I2OO2I : iLiebert et al.l l2003l : IWest et al. 



2004 ). More interestingly, we clearly find that the level of variability increases with later 
spectral type, or equivalently with decreasing log(LHa/ivboi)- This result is similarly evident 
from a comparison of the ratio of maximum to minimum Ha luminosity to the mean Ha 
luminosity (Figure [7]). 
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4. Individual Variability Events and Timescales 

To study the distribution of variability amplitudes and timescales in greater detail, we 
identify all of the individual variability "events" from the light curves; we do not include 
the non-variable light curves in this analysis. Events are defined as EW peaks that rise by 
at least 3o" above the nearest troughs. For example, in Figure [3fe, we identify three events 
— a broad event that lasts from about 10 to 70 min, and two events that last < 10 min 
at 20 and 35 min. Similarly, we find four events in Figure [3)d — a decline in the first 15 
min, followed by two short spikes, and finally a gradual rise between 40 and 65 min. The 
timescale of each event is defined as the time for the Ha EW to transition from a trough 
through a statistically significant peak to the next trough. In the case of partial events we 
use the observed timescale and amphtudes as lower limits. We find 71 full events and 27 
partial events in the light curves shown in Figures [1] and [2l 

A histogram of the variability event timescales is shown in Figure [S^. From the raw 
event list we find a peak at about 30 min, while partial events are naturally clustered at 
shorter timescales. The relatively small number of events with durations of ~ 10 min reflects 
a real trend since such events can be easily detected in our light curves. On the other hand, 
the decline beyond 30 min, corresponding to timescales longer than about one half of an 
observing sequence, may reflect the diminishing probability of capturing full events with a 
duration similar to that of the observation. For example, for our typical observations (1 hr 
duration with a time resolution of 5 min) the probability of detecting a full 40 min event is 
only 5/9 of the probability of detecting a full 20 min event. To take this effect into account 
we normalize each full event bin by its relative probability (Figure [8]d). 

In addition, to account for the distribution of partial events we make the simple as- 
sumption that these events have a uniform probability across all bins with a duration equal 
to or greater than their measured duration. This assumption is somewhat simplistic given 
the observed non-uniform trend in the full events, and it is thus likely to under-estimate the 
true number of events in the 30 — 50 min bins. However, given the overall number of events 
this is unlikely to change our conclusions in a significant way 

We show the histogram corrected for event detection probability and taking into account 
the partial events in Figure [8]d. As expected, the trend of increased frequency of events 
between timescales of 10 to 30 min becomes significantly more pronounced, with a factor of 
6 times as many events with 30 min duration compared to 10 min duration. The decline at 
longer timescales is shallower compared to the raw data, with a decline of about 30% relative 
to the 30 min bin (compared to 60 — 90% in the raw distribution). Taking into account our 
simple prescription of assigning partial events, it is likely that the distribution is in fact flat 
on timescales of 30 — 60 min. 
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We next turn to the event amplitude ratios. We focus on this quantity since it can 
be uniformly compared across spectral types. We plot the histogram of i?(EW) values in 
Figure M The distribution exhibits an exponential decline in the number of events (A^) as a 
function of amplitude ratio, with 

iVc.exp(-^2^). ,1) 

We note that this fit does not include partial events. To investigate whether there is a 
difference in the distribution as a function of spectral type, we divide the sample into two 
subsets split at M6 and repeat the analysis of timescale and amplitude distributions. We 
find the two sub-samples to be indistinguishable, indicating that there no obvious trend in 
Ha variability amplitude and timescale for individual events. 

Finally, we search for a correlation between the event durations and amplitudes. Fig- 
ure [10] shows the event amplitude ratios plotted against their durations. We find that R{EW) 
generally increases with increased timescale, but this effect is only apparent on timescales 
shorter than about 30 min. Longer events tend to have lower amplitudes. However, partial 
events are likely to increase the number of long duration, high amplitude events, potentially 
leading to a flatter distribution beyond 30 min. 



5. Discussion and Conclusions 

We carried out spectroscopic observations of 43 M dwarfs in the range M3.5-M8.5, 
for approximately 1 hour each with a time resolutions of about 5 min. About 80% of our 
targets exhibit statistically significant Ha variability, ranging from a factor of 1.25 to about 
4. Based on a total on-source exposure time of about 3000 min we find that the duty cycle 
of flares with an order of magnitude increase in brightness is < 0.05 hr~^ (95% confidence 
level, assuming Poisson distribution ). This limit is similar to previous results for individual 



well-studied objects (e.g., LHS2065: Martin fc Ardilall200lh . 



The level of variability for individual objects is found to increase with later spectral 
type, with a possible flattening beyond ~ M7. In terms of Ha luminosity, we recover the 
familiar trend of decreasing Lno/i^boi with later spectral type. More importantly, however, 
we find that the range of light curve variability increases with later spectral type, such that 
M7-M8 objects exhibit a range of about 0.5 dex in L^^a/Lhoi compared to only 0.15 dex for 
M4-M5 objects. This result indicates that while the traditional definition of "activity" (i.e., 
Lna/Lhoi) declines with later spectral type, the actual fluctuations in activity increase with 
later spectral type. 
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The individual light curves exhibit a rich phenomenology, with activity timescales span- 
ning the full range covered by our observations, i.e., ~ 5 — 60 min. We find that fluctuations 
on timescales of ~ 10 min are signiflcantly less common than those on ~ 30 min timescale, 
and that longer duration events (at least to ~ 1 hr) are likely to be as common. The 
event amplitude ratios closely follow an exponential distribution with a characteristic value 
of i?(EW) — 1 0.7. Taken at face value, this means that events with an amplitude ratio 
of 10 are ~ 10^ times less common than those with a ratio of 2. Thus, our observed rate 
of about 0.1 event per hour with a factor of 2 increase in EW implies an expected rate of 
~ 10^^ hr^^ for events with an order of ma gnitude increase in EW. The fact that several such 



events are published in the literature (e.g.. lMartm fc Ardilall200ll ) suggests that major flare 



events are not drawn from the same statis tical distribution of mild variability events found 



in our work. Indeed, ISchmidt et al.l (120071 ) estimate that events with A(EW) > 15 A have 
a duty cycle of ~ 5% in late-M dwarfs, orders of magnitude larger than expected from our 
exponential distribution (however, A(EW) depends on spectra type so a direct comparison 
is challenging). 

Combining the event amplitudes and timescales, we do not flnd any clear correlations, 
although there appears to be a general trend of declining amplitude for longer duration 
events. Such a relation would be expected if Ha events in each spectral type released similar 
amounts of total magnetic energy such that E'eq ~ x t oc EW x t ~ const. On the other 
hand, the opposite trend would be expected if Ha events represent the release of magnetic 
stresses in the chromosphere such that frequent and/or shorter duration events will have 
lower amplitude. It remains to be seen from observations with longer time baselines which 
effect exists, and whether it correlates with spectral type. 

The large and growing samples of nearby M dwarfs are conducive for a continued system- 
atic investigation of Ha emission as a proxy for magnetic activity. Clearly, Ha variability at 
the level of < 2 is prevalent in the bulk of Ha-emitting mid- and late-M dwarfs on timescales 
of ~ 0.1 — 1 hr. Future observations will need to address three primary questions: 



1. Are the durations and amplitudes of Ha events correlated? If so, directly or inversely? 

2. Are large flares drawn from the same distribution as small events? 

3. Are the prevalence and properties of events/flares correlated with age or rotation ve- 
locity across the M spectral type range? 



To address these questions we will continue to pursue observations of larger samples, as well 
as longer time baseline observations of individual active objects. 
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Table 1. Source Properties and Log of Observations 



Source 




Sp. Type 


Date=- 
(UT) 


Exposures log(Lboi/L0 


) " 


V sin i 
^^Kin S ) 


xxei. 


G 99-049 


GJ 3379 


M3.5 


Mar 


15 


23:58:23 


12x300 


-1.63 


5.4 


7.4 


1, 8 


LHS 1723 




M4 


Mar 


14 


00:01:34 


11x300 


-1.92 


6.1 


<3.2 


2, 8 


L 449-1 




M4 


Mar 


14 


01:00:14 


11x300 




5.7 




3 


GJ 1224 




M4.5 


Mar 


14 


08:29:12 


14x300 


-2.36 


7.5 


<5.6 


2, 8 


GL 285 


V* YZ CMi 


M4.5 


Mar 


15 


01:10:51 


13x300 


-1.52 


6.2 


6.5 


2, 8 


2MASSW J1013426-275958 




M5 


Maj: 


14 


04:28:58 


6x600 








4, 8 


GJ 1156 


V* GL Vir 


M5 


Mar 


14 


05:57:50 


12x300 


-2.30 


6.5 


9.2 


2, 8 


GJ 1154A 




M5 


Mar 


15 


03:41:37 


18x300 




8.5 


5.2 


2, 8 


DENIS-P J213422. 2-431610 




M5.5 


Sep 


12 


04:12:02 


4x300 




14.6 




5 


2MASS J02591181+0046468 




M5.5 


Sep 


14 


04:48:39 


11x300 


-3.61 


29 




6 


2MASS J02534448- 7959133 




M5.5 


Sep 


14 


06:21:38 


13x300 


-3.44 


17.2 




7 


2MASS J00244419-2708242 


GJ 2005 


M5.5 


Sep 


14 


07:42:30 


12x300 


-2.59 


7.5 


9.0 


8, 8 








Sep 


15 


03:20:25 


10x300 
















Sep 


15 


06:14:05 


38 X 300 










2MASS J00045753-1709369 




M5.5 


Sep 


15 


00:52:58 


9x300 


-3.31 


14.9 




5 


2MASS J20021341-5425558 




M5.5 


Sep 


12 


02:51:06 


11x300 


-3.55 


17.2 




7 


LP 844.25 


LHS 2067 


M6 


Mar 


14 


03:19:25 


6x600 


-3.90 


25.1 




9 


2MASS J16142520-0251009 


LP 624-54 


M6 


Maj: 


17 


08:34:42 


15x300 


-3.42 


14.6 




7 








Sep 


13 


23:47:29 


10x300 










2MASS J21322975-0511585 


NLTT 51488 


M6 


Sep 


14 


23:52:34 


6x300 


-3.47 


18.5 




7 


2MASS J23373831-1250277 


NLTT 57439 


M6 


Sep 


14 


01:47:33 


16x300 


-3.49 


19.2 




10 


2MASSW J1012065-304926 




M6 


Mar 


17 


01:13:30 


6x600 








4 


LP 731-47 




M6 


Mar 


17 


04:29:30 


12x300 




20.5 


11.0 


7, 8 


2MASS J23155449-0627462 


NLTT 56283 


M6 


Sep 


16 


4/\ Af\ 

00:40:42 


9x300 


-3.37 


17.7 




11 


2MASS J20424514.0500193 


NLTT 49734 


M6.5 


Sep 


13 


03:32:18 


12x300 


-3.51 


15.5 




12 


GJ 3622 




M6.5 


Mar 


15 


05:25:04 


12x300 




4.5 


3.0 


1, 8 


2MASS J05023867-3227500 




M6.5 


Mar 


16 


23:56:31 


12x300 


-3.89 


25.1 




10 


2MASS J02141251-0357434 


LHS 1363 


M6.5 


Sep 


13 


09:01:56 


10x300 


-3.11 


10.1 




10 








Sep 


14 


03:25:22 


4x300 










2MASS J10031918-0105079 


LHS 5165 


M7 


Mar 


14 


02:08:39 


6x600 


-3.84 


23.1 




10 


2MASS J13092185-2330350 




M7 


Mar 


14 


07:11:58 


6x600 


-3.62 


13.3 


7.0 


4, 13 


2MASSW J1032136-420856 




M7 


Mar 


15 


02:42:22 


6x600 








4 


2MASSW J1420544-361322 




M7 


Mar 


15 


06:37:15 


8x300 








4 


2MASS J09522188-1924319 




M7.5 


Mar 


16 


02:26:22 


10x450 


-3.68 




6.0 


13, 13 








Mar 


17 


23:45:15 


8x600 










2MASS J04291842-3123568 




M7.5 


Sep 


12 


07:12:29 


12x300 


-3.26 


11.4 




14 


2MASS J23062928-0502285 




M7.5 


Sep 


13 


05:11:27 


6x600 


-3.46 


11.0 




10 


2MASS J03313025-3042383 


NLTT 11163 


M7.5 


Sep 


13 


06:28:07 


6x600 


-3.46 


12.1 




10 


2MASS J04351612-1606574 


NLTT 13580 


M7.5 


Sep 


13 


07:44:12 


12x300 


-3.09 


8.6 




15 


2MASS J06572547-4019134 




M7.5 


Mar 


16 


01:22:06 


6x600 


-4.01 


22.7 




10 


2MASS J05173766-3349027 




M8 


Sep 


12 


08:27:48 


9x600 


-3.71 


14.7 




10 


2MASS J19165762-I-0509021 


VB 10 


M8 


Sep 


13 


00:54:33 


12x300 


-2.88 


5.7 


6.5 


8, 8 


2MASS J22062280-2047058 




M8 


Sep 


14 


00:35:30 


6x600 


-3.88 


18.2 


22.0 


10, 8 








Sep 


15 


23:39:08 


6x600 










2MASS J02484100-1651216 




M8 


Sep 


14 


09:10:55 


4x600 


-3.94 


16.2 




10 



Table 1 — Continued 



Source 


Other Sp. Type 




Date'' 


Exposures 


log(ibol/-f'0) 


b d 


V sin i Ref . 








(UT) 


(s) 




(pc) 


(km s^-"-) 


2MASS J20370715-1137569 


M8 


Sep 


16 23:39:45 


5x600 


-3.S5 


16.8 


10 


2MASS J22264440- 7503425 


M8.5 


Sep 


12 01:26:16 


6x600 


-3.S7 


16.5 


5 






Sep 


15 04:35:01 


6x600 








2MASS J03061159-3647528 


MS. 5 


Sep 


12 04:50:23 


5x600 


-3.61 


11.3 


5 


2MASS J23312174-2749500 


MS. 5 


Sep 


13 02:13:20 


6x600 


-3.61 


11.6 


5 



^AU observations were carried out in the year 2007. 

''Bolometric luminosities were derived usi ng bolometric correc tions on the J and K magnitudes listed in the 
SiMBAD database, using the fits described in IWilking et ail (Il999l') . 



'^Reference s for 
iHenrv et aT~l Il99 



spect ral types are i n stan dard font: 



■ ■ - - l) : [21 iDelfosse et al.l lll998h: 
[61 iBochanski et al.l 120051) ; [7] | Phan-Bao fcBessell 



J31 



Scholz et al 
200^ 



referen ces fo r rotational 
2005ah: 



veloc ities are in ita lics: 

■ _ _ _ ■20021 '): [ 51 ICrifo et al, 

[81 iMoharitv fc Basril mO±-. [91 Ifieid fc Gizi; 



[41 iGizii 



[lOllCruz et al.lll200ai): [UllScholz et alj (|2005bl) ril2] iReid et al.l ||2004| ') : [13] lReid et^ hom ; [14] lschmidt et all 
ll2007h : [15l lLodieu et all l l2005h . 



11] 



12005): 
1 2OO5I') 



Table 2. Ha Variability 



Source Sp. Type Published EW Ref.'' (EW)'' Min(EW) Max(EW) A(EW)'^ RMS(EW)'* Mean Max 

log(LHc«/ibol) logCiHa/i'bol) 







(A) 




(A) 


(A) 


(A) 


(A) 


(A) 






G 99-049 


M3.5 


2.9 


1 


7.1 ± 0.1 


7.0 ± 0.1 


8.1 ± 0.1 


1.1 ± 0.1 


0.5 ± 0.1 


-3.229 


-3.178 


LHS 1723 


M4 


0.9 


1 


1.2 ± 0.2 


1.0 ± 0.1 


1.3 ± 0.1 


0.3 ± 0.2 


0.1 ± 0.1(C) 


-4.155 


-4.082 


L 449-1 


M4 






6.8 ± 0.2 


6.7 ± 0.1 


7.1 ± 0.1 


0.4 ± 0.2 


0.1 ± 0.1(C) 


-3.357 


-3.345 


GJ 1224 


M4.5 


2.3 


1 


4.8 ± 0.2 


4.3 ± 0.2 


5.2 ± 0.2 


0.9 ± 0.2 


0.2 ± 0.2 


-3.480 


-3.436 


GL 285 


M4.5 


9.5 


1 


9.3 ± 0.1 


8.6 ± 0.1 


10.6± 0.1 


2.0 ± 0.1 


0.6 ± 0.1 


-3.179 


-3.152 


2MASSW J1013426-275958 


M5 


6.6 


2 


7.2 ± 0.4 


6.5 ± 0.4 


7.8 ± 0.3 


1.3 ± 0.5 


0.4 ± 0.3(C) 


-3.711 


-3.670 


GJ 1156 


M5 


4.4 


1 


6.9 ± 0.3 


6.4 ± 0.2 


8.8 ± 0.2 


2.3 ± 0.3 


0.7 ± 0.2 


-3.705 


-3.618 


GJ 1154A 


M5 


4.3 


1 


8.9 ± 0.3 


6.8 ± 0.2 


13.5± 0.2 


6.6 ± 0.2 


2.0 ± 0.2 


-3.589 


-3.448 


DENIS-P J213422.2-431610 


M5.5 






0.9 ± 0.5 


0.6 ± 0.3 


1.6 ± 0.3 


1.1 ± 0.4 


0.4 ± 0.4(C) 


-4.466 


-4.262 


2MASS J02591 181+0046468 


M5.5 


14.5 


3 


12.0 ± 0.8 


9.3 ± 0.7 


15.2± 0.4 


5.8 ± 0.8 


1.9 ± 0.6 


-3.387 


-3.290 


2MASS J02534448- 7959133 


M5.5 


12.4 


4 


7.8 ± 0.3 


6.3 ± 0.4 


12.2± 0.6 


5.9 ± 0.7 


1.9 ± 0.5 


-3.502 


-3.387 


2MASS J00244419-2708242 


M5.5 


3.6 


1 


4.0 ± 4.5 


1.8 ± 4.5 


8.5 ± 0.3 


6.7 ± 4.5 


1.1 ± 0.7 


-3.734 


-3.537 


2MASS J00045753-1709369 


M5.5 






4.4 ± 0.4 


3.0 ± 0.3 


5.9 ± 0.4 


2.9 ± 0.5 


0.8 ± 0.4 


-3.813 


-3.695 


2MASS J20021341-5425558 


M5.5 


7.6 


4 


1.9 ± 0.4 


3.5 ± 0.8 


7.5 ± 0.5 


4.0 ± 0.9 


1.3 ± 0.6 


-3.785 


-3.591 


LP 844-25 


M6 






6.7 ± 1.2 


1.0 ± 0.3 


2.3 ± 0.3 


1.2 ± 0.4 


0.5 ± 0.3 


-4.610 


-4.394 


2MASS J16142520-0251009 


M6 


4.2 


4 


8.5 ± 0.9 


1.9 ± 0.7 


5.9 ± 0.6 


3.9 ± 0.9 


1.0 ± 0.6 


-4.187 


-3.985 


2MASS J21322975-0511585 


M6 


1.1 


4 


3.5 ± 0.8 


0.1 ± 0.5 


1.4 ± 0.5 


1.4 ± 0.7 


0.6 ± 0.5 


-5.057 


-4.610 


2MASS J23373831-1250277 


M6 






3.1 ± 0.7 


14.6± 0.5 


25. 9± 0.6 


11. 3± 0.8 


4.9 ± 0.5 


-3.500 


-3.358 


2MASSW J1012065-304926 


M6 


5.6 


2 


12.5 ± 0.7 


4.6 ± 0.7 


10.4± 0.4 


5.8 ± 0.8 


2.2 ± 0.7 


-3.904 


-3.756 


LP 731-47 


M6 


6.4 


4 


0.6 ± 0.7 


4.8 ± 0.7 


11. 4± 1.0 


6.6 ± 1.2 


2.0 ± 0.7 


-3.847 


-3.714 


2MASS J23155449-0627462 


M6 


1.1 


5 


5.0 ± 0.4 


3.4 ± 0.9 


5.9 ± 0.7 


2.5 ± 1.2 


1.1 ± 0.8 


-4.122 


-3.985 


2MASS J20424514^0500193 


M6.5 






4.1 ± 0.4 


2.1 ± 0.5 


4.5 ± 0.4 


2.4 ± 0.6 


0.6 ± 0.5 


-4.451 


-4.341 


GJ 3622 


M6.5 


1.9 


1 


7.7 ± 0.8 


3.4 ± 0.3 


4.6 ± 0.3 


1.1 ± 0.4 


0.3 ± 0.3(C) 


-4.393 


-4.332 


2MASS J05023867-3227500 


M6.5 






3.3 ± 0.7 


7.4 ± 0.5 


12.3± 2.5 


4.9 ± 2.5 


2.5 ± 1.0 


-4.050 


-3.915 


2MASS J02141251-0357434 


M6.5 






7.2 ± 0.4 


6.9 ± 0.3 


18. 3± 0.5 


11. 5± 0.6 


3.7 ± 0.4 


-4.012 


-3.732 


2MASS J10031918-0105079 


M7 






10.6 ± 1.2 


8.9 ± 0.9 


17.4± 0.7 


8.4 ± 1.1 


2.9 ± 0.8 


-4.201 


-4.050 


2MASS J13092185-2330350 


M7 


6.7 


2 


6.0 ± 0.9 


5.1 ± 0.6 


9.3 ± 0.5 


4.2 ± 0.8 


1.5 ± 0.6 


-4.383 


-4.312 


2MASSW J 1032 136-420856 


M7 


11.8 


2 


11.9 ± 1.3 


11.9± 1.3 


19.9± 1.4 


8.0 ± 1.9 


3.7 ± 1.4 


-4.201 


-4.002 


2MASSW J1420544^361322 


M7 


7.0 


2 


16.5 ± 1.2 


8.5 ± 0.6 


30.0± 0.5 


21.6± 0.8 


8.5 ± 0.8 


-4.050 


-3.804 


2MASS J09522188-1924319 


M7.5 


11.3 


6 


4.5 ± 2.1 


8.1 ± 0.6 


18.9± 0.3 


10.8± 0.7 


3.1 ± 0.4 


-3.940 


-3.736 


2MASS J04291842-3123568 


M7.5 


15.9 


7 


10.9 ± 0.7 


10.8± 0.4 


14.8± 0.5 


3.9 ± 0.6 


1.2 ± 0.5 


-3.933 


-3.866 


2MASS J23062928-0502285 


M7.5 


2.8 


7 


12.6 ± 0.5 


3.4 ± 0.7 


5.4 ± 0.8 


2.0 ± 1.1 


0.8 ± 0.7(C) 


-4.379 


-4.280 


2MASS J03313025-3042383 


M7.5 


7.6 


7 


3.5 ± 0.9 


7.4 ± 0.6 


10.9± 0.6 


3.5 ± 0.8 


1.5 ± 0.7 


-4.068 


-4.012 



Table 2 — Continued 





Source 


bp. iype 


r^UDlisnca rLYv 


ixei. 


{^rj W / 


lViin(^£!j W ) 


IViaxi^rljW ) 


Zi^rL/W ) 


KiViol^ilf W ) 


iviean 


Max 
























Jogl^Hd/^bolj 








(A) 




(A) 


(A) 


(A) 


(A) 


(A) 






2MASS 


J04351612-1606574 


M7.5 


3.7 


6 


7.8 ± 1.0 


5.3 ± 0.8 


8.3 ± 0.4 


2.9 ± 0.8 


0.9 ± 0.6 


-4.206 


-4.093 


2MASS 


J06572547-4019134 


M7.5 






5.7 ± 0.8 


4.4 ± 1.2 


11. 8± 1.7 


7.4 ± 2.1 


2.5 ± 1.8 


-4.180 


-3.971 


2MASS 


J05173766-3349027 


M8 


8.1 


4 


2.7 ± 1.6 


1.5 ± 1.2 


5.4 ± 1.1 


3.9 ± 1.6 


1.3 ± 1.1(C) 


-4.688 


-4.487 


2MASS 


J19165762+0509021 


M8 


5.6 


1 


4.7 ± 0.8 


3.8 ± 0.5 


8.8 ± 0.5 


5.1 ± 0.7 


1.3 ± 0.5 


-4.479 


-4.275 


2MASS 


J22062280-2047058 


M8 


5.1 


8 


3.6 ± 1.3 


2.3 ± 0.9 


12.5± 1.2 


10.2± 1.5 


2.9 ± 1.1 


-4.356 


-4.123 


2MASS 


J02484100-1651216 


M8 


7.9 


7 


6.5 ± 1.2 


2.1 ± 1.3 


7.4 ± 1.0 


5.3 ± 1.7 


1.9 ± 1.2 


-4.538 


-4.350 


2MASS 


J20370715-1137569 


M8 


6.2 


7 


1.8 ± 1.3 


0.7 ± 0.7 


2.1 ± 1.4 


1.4 ± 1.5 


0.5 ± 1.2(C) 


-5.016 


-4.897 


2MASS 


J22264440-7503425 


M8.5 


3.2 


7 


0.6 ± 2.4 


3.3 ± 1.1 


11. 7± 1.5 


8.4 ± 1.9 


2.7 ± 1.7(C) 


-4.538 


-4.350 


2MASS 


J03061159-3647528 


M8.5 


0.5 


6 


8.3 ± 0.7 


3.9 ± 0.8 


12.1± 0.7 


8.1 ± 1.1 


2.6 ± 0.8 


-4.421 


-4.239 


2MASS 


J23312174-2749500 


M8.5 


5.1 


6 


5.5 ± 1.0 


3.5 ± 0.9 


7.1 ± 0.7 


3.6 ± 1.1 


1.1 ± 0.9 


-4.594 


-4.467 



^"Re fere nces for the publishe d E Ws: fDlMohantv fc Basril ll2003h f2')lGizij ||2002| ') f iBochanski et al.ll|2005h f4') |Phan-Bao fc Besselil l|2006l ') f5') lReid et al] 
l|2003tl (6) iLodieu et all tood ) (7) ISchmidt et al.l fcoOTl' l (Sl lReid et al.l ll2002li 

•"Median EW. 

<=A(EW) = Max(EW) - Min(EW) 

"^'(C)' denotes that the object's EW was constant within errors during the observations (see discussion in SectionlSj. Note that while 2M2226-7503 had 
two separate observations, both lightcurves were found to be consistent with a constant EW. 
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G099-049, M3.5 
I ' ' ' I 



LHS1723, M4.0 



U49-1, M4.0 



GJ1224. M4.5 



GL285. M4.5 




20 40 60 
Mins after MJD54 174.99887 



20 40 
Mins after MJD54174.00I06 



10 20 30 40 
Mins after MJDS4174.0418S 



20 40 60 80 
HIns after MJD54174.35361 



20 40 60 
Mins after MJD54 175.04920 



2M1013-Z759, M5.0 



GJ1156, M5.0 



GJ1154A, M5.0 



2M2134-4316, M5.5 



2M0259 + 0046, M5.5 

I- I I I I I I I I I I 




20 40 
Mins after MJDS4 174. 16678 



20 40 60 
Mins after MJDS4 174.34849 



50 100 
Mins after MJD54175.15390 



5 10 15 
Mins after MJD54356. 17502 



20 40 60 
Mins after MJD543e8.2004& 



2M0253-7959. M5.5 



2M0004-1709. M5.5 

TJTTTTJTTT 



2M0024-2708, M5.5 



2M0024-2708, M5.5 



2M0024-2708. M5.5 




20 40 60 10 20 30 40 20 40 60 20 40 100 200 

Mins after MJD543S8.26502 Mins after MJD543&9.Q3678 Mins after MJD543S8.32UB Mins after MJD&43S9.13917 Mins after MJD54359.2597B 



LP844-25, M6.0 



2M1012-3049, M6.0 



LP731-47, M6.0 



2M1614-0251, M6.0 



2M1614-0251, M6.0 




20 40 



20 40 60 



20 40 60 



SO 40 60 80 



20 40 



Mins after HJD54174.I384B Mins after UJDS4177.05104 Uins after MJDS4177.I8715 Mins after IUD54 177.35743 Mins after HJD5435e.99130 



2M2002-5426, M6.0 




25 



•< 20 



15 



2M2337-1250, M6.0 

.1 1 1 1 1 1 1 ijj 11 1 1 1 1 1 




2 
1.5 
1 

0.5 


-0.5 



2M2132-0511, M6.0 




20 40 60 



20 40 60 80 



i_Li 



2M2315-0627, M6.0 



GJ3622, M6.5 



10 20 30 




10 20 30 40 



20 40 60 



Uins after HJD543Se.l 1881 Mins after 1(4054358.07468 Mins after HJDS4357.99483 Mins after UJD543S0.028S6 Mins after MJDS4175.SSS74 



Fig. 1. — Ha EWs plotted against time in minutes, for all observations in our sample. Note 
that some objects have several light curves from repeated observations. 
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Mins after MJD54 175.99758 Mins after MJD54357. 14743 Mins after MJD543S7. 37634 Mlns after MJD54358. 14261 Mins after MJD54174.0B934 




20 40 10 20 30 40 10 20 30 40 20 40 60 BO 20 40 60 80 

Mins after MJD54 174 .29997 Mins after MJD54175.11275 Mins after MJD54175.2758e Mins after MJD54 176. 10 164 Mins after MJD54 176.98975 



3M0657-4019, M7.5 3M0429-3134, M7.5 aM2306-0503, M7.5 3M033 1-3043, M7.5 3M0435-1606, M7.5 




30 40 60 60 30 40 60 30 40 30 40 10 30 30 

Mins after MJDS43Se.35S63 Mins after MJD543&7.D37e8 Mins after MJD54358.02465 Mins after MJD54358.9B&50 Mins after MJD543S8.36258 



2M2037-1137. M6.0 2M2226-7503. M8.5 2M2326-7503. M6.5 2M0306-3647. M8.5 2M2331-2749, M8.5 




10 20 30 40 20 40 20 40 50 100 20 40 

Mins after MJD54359.9e593 Mins after MJD54356.05990 Mins after MJD5435g.l9098 Mins after MJD54356.30ieS Mins after MJD54357.09Z59 



Fig. 2. — Same as Figure [H 
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Fig. 3. — Subset of Figures [T] and [2] enlarged for discussion. Figure [3^ is an example of a 
light curve consistent with non- varying EW (see Figures [3]d and [31:; display examples 
of objects with rapid and gradual variability, respectively. The other plots show 3 of the 
objects that displayed the greatest variability during our observations. 
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Fig. 4. — Left: Difference between the maximum and minimum Ha EW, A(EW), plotted 
against spectral type. Center: Normalized RMS of the Ha EW light curves of individual 
observations, plotted against spectral type. Right: Ratio of maximum and minimum Ha 
EW, R(EW), plotted against spectral type. In all three plots the error bar on the upper left 
show the median errors. Open circles designate objects identified as non-varying using our 
criterion (Q. 
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5 10 15 

<EW> (A) 



Fig. 5. — Median-normalized RMS of the EW light curves plotted against the median EW. 
The error bars at upper left show the median errors. Open circles denote objects identified 
as non- varying through our criterion. 
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Fig. 6. — Logarithmic ratios of and Lboi, plotted against spectral type. The solid 
lines connect the maximum and minimum i^Ha/^boi values measured for each object. Open 
circles designate objects identified as non- varying through our criterion. The positions 
of objects arc displaced horizontally for clarity, and the dashed vertical lines delineate the 
spectral types. 
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Fig. 7. — Ratio of maximum and minimum Ha luminosity plotted against median Ha lumi- 
nosity. The error bar on the upper right show the median uncertainty. Open circles denote 
objects identified as non-varying through our criterion. 
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Fig. 8. — Left: Histogram of variability events binned by timescale. Unshaded regions denote 
partially observed events. Right: Frequency of variability events weighted by the probability 
of being observed within a finite observing time of about 1 hr, and with the partial events 
uniformly distributed on timescales equal to or longer than their observed timescale. The 
normalization is chosen such that the total area is unity. 
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Flare amplitude ratio: Max(EW) / Min(EW) 



Fig. 9. — Histogram of variability events binned by flare amplitude ratio. The solid curve is 
an exponential distribution with a decay constant of 0.7. Unshaded regions denote partially 
observed events. 
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Fig. 10. — Amplitude ratio of variability events plotted against timescale. Square boxes with 
diagonal arrows denote partially observed events, which are lower limits in both amplitude 
and timescale. The error bars at the upper right show the median errors. 



